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in optoelectronic oscillators. After a short overview of the main results 
for the phase diffusion approach in autonomous oscillators, an extension 
is proposed for the case of optoelectronic oscillators where the microwave 
is a limit-cycle originated from a bifurcation induced by nonlinearity and 
time-delay. This Langevin approach based on stochastic calculus is also 
successfully confronted with experimental measurements. 
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1 Introduction 

Optoelectronic oscillators (OEOs) combine a nonlinear modulation of laser light 
with optical storage to generate ultra-pure microwaves for lightwave telecom- 
munication and radar applications [1, 2]. Their principal specificity is their 
extremely low phase noise, which can be as low as —160 dBrad 2 /Hz at 10 kHz 
from a 10 GHz carrier. Despite some interesting preliminary investigations, the 
theoretical determination of phase noise in OEOs is still a partially unsolved 
problem. The qualitative features of this phase noise spectrum can be recov- 
ered using some hcuristical guidelines or rough approximations, but however, a 
rigorous theoretical background is still lacking. 

There are several reasons which can explain that absence of theoretical back- 
ground. A first reason is that before refs. [3], there was no time domain model 
to describe such systems, so that stochastic analysis could not be used to per- 
form the phase noise study. Moreover, unlike most of oscillators, the OEO is a 
delay- line oscillator, and very few had been done to study the effect phase noise 
on time-delay induced limit-cycles. Finally, the OEO is subjected to multiple 
noise sources, which are sometimes non-white, like the flicker (also referred to 
as "1//") noise which is predominant around the microwave carrier. 

The objective of this work is to propose a theoretical study where all these 
features are taken into account. The plan of the article is the following. In 
Section 2, we present the phase diffusion approach in autonomous systems. It 
is a brief review where the fundamental concepts of phase diffusion are recalled, 
and where some important earlier contributions are highlighted. Then, we derive 
in Section 3 a stochastic delay-differential equation for the phase noise study. 
We show that for our purpose, the global interaction of noise with the system can 
be decomposed into two contributions, namely an additive and a multiplicative 
noise contribution. Section 4 is devoted to the study of the noise spectrum below 
threshold. It will appear that the spectrum below threshold will not only be 
important to validate the stochastic model, but also that it enables an accurate 
calibration of additive noise. In Section 5, we address the problem of phase 
noise when there is a microwave output using Fourier analysis, and we show 
that it is possible to have an accurate image of the phase noise spectrum in all 
frequency ranges. The last section concludes the article. 

2 The phase diffusion approach in autonomous 
oscillators 

2.1 Fundamental concepts 

For an ideal (noise-free) oscillator, the Fourier spectrum is a collection of Dirac 
peaks, standing for the fundamental frequency and its harmonics. The effect of 
amplitude white noise is to add a flat background, while the peaks do keep their 
zero linewidth; it is the effect of phase noise to widen the linewidth of these 
peaks. 

Some pioneering papers on the topic of phase noise in autonomous oscillators 
using stochastic calculus had been published forty years ago [4]. In particular, it 
was demonstrated that a general framework to study the problem of phase noise 
in a self-sustained oscillator could be built using some minimalist assumptions. 
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Figure 1: Experimental set-up. 



The first point is that a strong nonlinearity is an essential necessity in oscillators, 
in the sense that nonlinearity can not be regarded as small because it controls 
the operating level of the oscillator. The second important point is that the 
phase is only neutrally stable, so that quasilinear methods which assume that 
fluctuations from some operating point are small (linearization techniques) can 
not be applied directly. 

The phase is neutrally stable as a consequence of the phase-invariance of 
autonomous oscillators. In other words, limit-cycles are stable against amplitude 
perturbations, while there is no mechanism able to stabilize the phase to a 
given value: hence, phase perturbations are undamped, but they do not diverge 
exponentially, though. In a noise free oscillator, the "stroboscopic" state point 
on the limit-cycle is immobile, but in the presence of noise, it moves randomly 
along the limit-cycle: in other words, the phase of the oscillator undergoes a 
diffusion process, in all points similar to a one-dimensional Brownian motion. 
In the most simple case, the random fluctuations of the phase are referred to 
as a Wiener process, obeying an equation of the kind ip — where £ is a 
Gaussian white noise with autocorrelation (£(t)£(i')) — 2D5(t — t'), while D is 
a parameter referred to as the diffusion constant. It can be demonstrated that 
the phase variance diverges linearly as (<y9 2 (i)) — 2Dt, and the single-side band 
phase noise spectrum (in dBc/Hz) explicitly reads C(u>) = 2D/[D 2 + uj 2 ], so 
that D is the unique parameter characterizing all the statistical and spectral 
features of phase fluctuations. 

2.2 The unifying theory of Demir, Mehrota and Roychowd- 
hury 

On the base of earlier works by Lax [4] and Kartner [5], Demir, Mehrota and 
Roychowdhury have proposed few years ago a unifying theory of phase noise in 
self-sustained oscillators subjected to white noise sources [6]. Their approach, 
which had later been extended by Demir to the case colored noise sources [7], 
relies on stochastic calculus. The principal point of their contribution was the 
introduction of a decomposition of phase and amplitude noise through a pro- 
jection onto the periodic time-varying eigenvectors (the so called Floquet eigen- 
vectors; also see ref. [8]), and they proved that it provides the correct solution 
to the problem. 

Demir et al have shown that if the sources of noise are Gaussian and white, 
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the phase noise around the fundamental peak (and its harmonics) has a Lorentzian 
lineshape, and therefore is fully determined by an "effective" diffusion constant 
D e ff which is the unique parameter needed for the phase noise determination. 
However, if the Demir et al theory has the great and essential advantage of math- 
ematical rigorousness, its principal drawback is that exactitude is obtained at 
the expense of simplicity: the calculation of D e ff is very complex, as it requires 
an accurate determination of all the time-varying eigenvectors related to the 
autonomous flow. In general this task can only be performed numerically us- 
ing quite complicated algorithms, and this lack of flexibility explains why this 
method is scarcely used in the phase noise studies available in the literature. 
The key challenge for the study of phase noise in OEO would be provide an 
accurate description of the phase noise spectrum, while avoiding the determina- 
tion of Floquet eigenvectors, which is an extremely complicated task in delayed 
system. 

3 Application of the phase diffusion approach to 
OEOs: stochastic delay-differential equations 

The OEO under study is organized in a single-loop architecture as depicted 
in Fig. 1. The oscillation loop consists of: (i) A wideband integrated optics 
LiNb03 Mach-Zchndcr (MZ) modulator, seeded by a continuous-wave semicon- 
ductor laser of optical power P; the modulator is characterized by a half-wave 
voltage V-k = 4 V. (ii) A thcrmalized 4 km fiber performing a time delay of 
T = 20 ji on the microwave signal carried by the optical beam; the correspond- 
ing free spectral range is 0^/2^ = 1/T = 50 kHz. {Hi) A fast photodiode with 
a conversion factor S. (iv) A narrow band microwave radio-frequency (RF) 
filter, of central frequency Fo — Q,q/2tt = 10 GHz, and —3 dB bandwidth of 
AF = Ail/2iT — 50 MHz; (v) A microwave amplifier with gain G. (vi) A vari- 
able attenuator, in order to scan the gain, (vii) All optical and electrical losses 
are gathered in a single attenuation factor k. 

The dynamics of the microwave oscillation can therefore be described in 
terms of the dimensionless variable x(t) = TrV(t)/2V- R whose dynamics obeys [3] 

dx 1 f * 

x + t— + -/ x(s)ds = (3cos 2 [x(t-T) + cj)}, (1) 
dt V J t 

where /3 — ttkSGP^V^ is the normalized loop gain, <fr = ttVb/^V^ is the Mach- 
Zchnder offset phase, while r = 1/Af2 and 6 = Af2/f2g are the characteristic 
timcscalc parameters of the bandpass filter. Since we are interested by single- 
mode microwave oscillations, the solution of Eq. (1) can be expressed under the 
form 

x (t) = l -A{t) e lUot +^A* (t) e- <n °* , (2) 

where Ait) = A(t) exp[iip(t)] is the slowly varying amplitude of the microwave 
x(t). We can significantly simplify the right-hand side term of Eq. (1) because 
the cosine of a sinusoidal function of frequency Q can be Fourier-expanded in 
harmonics of f^o- I n other words, since x(t) is nearly sinusoidal around Oo, then 
the Fourier spectrum of cos 2 [a;(t — T) + 4>] will be sharply distributed around the 
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harmonics of using the relationship cos 2 z = [1 + cos2z]/2 and the Jacobi- 
Angcr expansion 



+oo 

e lzcosa = 2™JnO)e mQ , (3) 

n— — oo 

where J n is the n-th order Bessel function of the first kind. Hence, since the filter 
of the feedback loop is narrowly resonant around fio, it can be demonstrated 
that discarding all the spectral components of the signal except the fundamental 
is an excellent approximation, so that Eq. (1) can be rewritten as 

i dx n% /■* , . , n . „ 1 
x+ Andt + ^J to x{s)ds = -^ 

xJ![2\A(t - T)\] cos[fi (t-T)+V(t-T)]. (4) 

In order to include noise effects in this equation, we will consider two main noise 
contributions in this system. 

The first contribution is an additive noise, corresponding to random envi- 
ronmental and internal fluctuations which are uncorrelated from the eventual 
existence of a microwave signal. The effect of this noise can be accounted for 
by addition as a Langevin forcing term, to be added in the right-hand side of 
Eq. (4). This additive noise can be assumed to be spectrally white, and since 
we are interested by its intensity around the carrier frequency Clo, it can be 
explicitly written as 

Ut) = lUt)e iQot +lc(t)e- inot , (5) 

where ( a (t) is a complex Gaussian white noise, whose correlation is (( a (t)Q (*')) = 
AD a 8{t — t'), so that the corresponding power density spectrum is |£ a (w)| 2 = 
2D a . 

The second contribution is a multiplicative noise due to a noisy loop gain. 
Effectively, the normalized gain parameter explicitly reads 

V B 



n . „ , 7T KSGP . 

7 = /?sm2</> = -— sin 



(6) 



If all the parameters of the system are noisy (i.e., we replace k by K+5n(t), S by 
S+SS(t), etc.), then the gain 7 may be replaced in Eq. (4) by "f+Sj(t), where the 
5j(t) is the overall gain fluctuation. We therefore introduce the dimensionless 
multiplicative noise 

V m (t) = ^ , (7) 
7 

which is in fact the relative gain fluctuation. In the OEO configuration, we 
have rj m (t) <C 1. This noise is in general spectrally complex, as it is the sum of 
noise contributions which are very different (noise from the photodetector, from 
the amplifier, etc.). In agreement with the usual noise spectrum of amplifiers 
and photodetectors, we will here consider that this multiplicative noise is flicker 
(i.e., varies as l/f) near the carrier, and white above a certain knee- value. We 
therefore assume the following empirical noise power density 



\f, m (v)\ 2 =2D n 



u> + Q L 



(8) 
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where n^ is the low corner frequency of the flicker noise, while is the high 
corner frequency. More precisely we consider that the noise is white below 0,^ 
and above fin, while it remains flicker in between. Typically we may consider 
Q.l/2-k < 1 Hz and Uh/2tt > 10 kHz, so that the flicker noise is extended over 
a frequency span of more than 4 orders of magnitude. 

To avoid the integral term of Eq. (4) which is complicated to manage analyt- 
ically, it is mathematically convenient to use the intermediate integral variable 



i(t) = / x(s) ds = - 
Jt z 



B(t) e iUot + -B*(t)e- iQot . 



(9) 



which is also nearly sinusoidal with a zero mean value. Using Eqs. (4), (5) and 
(7) , it can be shown that the slowly- varying amplitude B(t) obeys the stochastic 
equation 

{B + (An + 2ifl )B + in An B}e mot + c.c. 



-2Att 7 [l + r] m (t)] 



I e ifio(t-T) e ^ T+c ^ 



xJ 1 [2\B T + in Q B T \] + 2An 



CaWe inot + c.c. 



(10) 



where c.c. stands for the complex conjugate of the preceding term. We can 
assume \B\ <C Af2|Z?| and \B\ <§; Oo|<B|; the relationship x(t) — u(t) therefore 
gives A ~ in B, so that we can finally derive from Eq. (10) the following 
stochastic equation for the slowly varying envelope A(t) 



A = -^A + 2 lt ie li >[l + T lm (t)] J Ci [2\At\]A t 



(11) 



where Jci(x) = Ji(x)/x is the first order Bessel cardinal function of the first 
kind. The phase condition has been set to e~ tn ° T = —1, so that the dynamics 
of interest is restricted to the case 7 > 0. The key parameters of this equation 
are 



An/2 



v/i + (i/2Q) 2 



and -d = arctan 



1 



2Q 



(12) 



where Q = n n /An = 200 is the quality factor of the RF filter. Since Q ^> 1, 
we may simply consider that /i ~ An/2 and ~ 1/2Q. The complex term 
\xe % ® is a kind of "filter operator", which can be simply equated to the half- 
bandwidth An/2 when the Q-factor of the filter is sufficiently high, as it was 
done in ref. [3]. It is also noteworthy that in the complex amplitude equation 
(11), the initial multiplicative noise remains a real variable, while the additive 
noise becomes complex. 

We had recently shown, in agreement with the experiment, that the OEO 
has three fundamental regimes [3]. For 7 < 1, the system does not oscillate 
and the trivial fixed point is stable; for 1 < 7 < 2.31, the system sustains a 
pure microwave oscillation, with a constant amplitude and frequency; and at 
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last, for 7 > 2.31, the system enters into a regime where the amplitude of the 
microwave is unstable, and turns to be nonlinearly modulated. We can consider 
that this phenomenology is still correct as long as Q » 1. With the aid of 
the stochastic delay-differential equation ruling the dynamics of A, we may now 
derive analytically the power spectrum density of the oscillator, below and above 
threshold. However, it should be stressed that in all cases, stochastic variables 
should be manipulated with respect to the rules of stochastic calculus when an 
integral/differential transformation is applied to them. 



4 Noise power density spectrum below thresh- 
old (7 < 1) 

In general, no interest is paid to the study of the noise power density spectrum 
below threshold in OEOs. This lack of interest can be explained by the fact 
that there is no oscillation in this regime, and the system randomly fluctuates 
around the trivial equilibrium. However, as we will further see, this regime is 
particularly interesting because it enables to understand how the noise interacts 
with the system. 

From the stability theory of delay-differential equations with complex coeffi- 
cients, the deterministic solution of Eq. (11) below threshold is the trivial fixed 
point A = 0. After linearization around this solution, Eq. (11) can simply be 
rewritten as 

A = -ne™ A + we™ A T + ^ Q{t) , (13) 

where we have used Jci(0) = 1/2. This equation indicates that the multiplica- 
tive noise has no significative influence below threshold, because the product 
rj m AT is a second-order term. Therefore, the noise power below threshold is 
essentially determined by additive noise. 

Equation (13) is linear with constant coefficients: hence, the power density 
spectrum can directly be obtained as 

|iM| 2 = (14) 

live-™ + n(l - -ye-™*') | 2 

One can determine the total output power below threshold due to the white 
noise fluctuations in the system through the formula 

where R is the output impedance (in our case, R = 50 ft). The dimensionless 
power (|^4(t)| 2 ) can not be calculated analytically for 7^0: it can nevertheless 
be determined cither by numerical simulation of Eq. (13), or through a numerical 
computation of the integral ^ |^4(w)| 2 dui, where A(uj) is given by Eq. (14). 

However, in the open-loop configuration (7 = 0), the noisy output power 
can be analytically determined as 

AV 2 

P =-^AFD ai (16) 
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Figure 2: Variation of the RF noise output power P 1 as a function of the 
normalized gain, under threshold. The solid line is the theoretical prediction 
of Eq. (15) with D a = 9.8 x 1CP 16 rad 2 /Hz, and the symbols represent the 
experimentally measured data. The gain was varied through attenuation in the 
electric branch of the loop. 

through the use the Fourier integral, or using fundamental results from stochas- 
tic calculus since Eq. (13) degenerates to the well-known Orstein-Uhlenbeck 
equation. Therefore, knowing the bandwidth AF of the RF filter and the half- 
wave voltage V w of the MZ interferometer, an open-loop measurement of the 
output power can directly give an experimental a value for the white noise 
power density D a through Eq. (16). 

In our system, we have experimentally measured Pq = 20 nW (or —47 dBm), 
which corresponds to D a = 9.8 x 1CT 16 rad 2 /Hz. This value for the power can 
also be obtained by other means [see Appendix A]. The curve displaying the 
power variation as a function of the normalized gain under threshold is shown in 
Fig. 2, and there is an excellent agreement between the experimental data and 
our analytical formula of Eq. (15). It may be interesting to note that the noise 
power apparently diverges at 7 = 1. In fact, one should not forget that this 
result is obtained using Eq. (13), which is only valid for |*4| <C 1. When 7 — > 1, 
the amplitude of A increases and the higher order terms of the Bessel cardinal 
function are not negligible anymore, so that the Eqs. (13) and (14) are no more 
valid. Hence, divergence of the noise power is prevented by the nonlinear terms 
of Eq. (11) which become predominant in a very narrow range just below the 
threshold. A noteworthy study on this topic of noisy oscillators near threshold 
isref. [10]. 

It is also noteworthy that for 7 = 0, the noise spectrum follows the spectral 
shape of the RF filter. However, when 7 is increased (still below threshold), a 
first qualitative difference emerges, since the spectrum still follows the spectral 
shape of the filter, but its fine structure is composed by a collection of peaks 
which are the signature of microwave ring-cavity modes, as it can be seen in 
Figs. 3 and 4. 
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Figure 3: Theoretical power density spectrum |„4(/)| 2 of the microwave noise 
signal below threshold, with D a — 9.8 x 1CP 16 rad 2 /Hz and 7 = 0.8, using 
Eq. (14). The semi-logarithmic scale is adopted because it enables to see at 
the same time the fine structure of regularly spaced ring-cavity peaks, and the 
global variation shaped by the RF filter bandwidth. This spectrum is divided 
into two areas: a quasi-flat area within bandwidth, and a —20 dB/dec decrease 
outside the bandwidth. 




5 Phase noise spectrum above threshold (7 > 1) 

Above threshold, the amplitude of the microwave obeys the nonlinear algebraic 
equation Jc 1 [2|^4 |] = 1/(27). Linearizing Eq. (11) around this solution yields 
the following equation 

A = -ii^A + /ie w [1 + Vm(t)] At + ^ d ( a (t) . (17) 

We should now look for an equation for the phase ip in order to find its power 
density spectrum ^(w)! 2 . 

Using the Ito rules of stochastic calculus [see Appendix B], we derive the 
following time-domain equation for the phase dynamics 

^ = -^-« + ^ m (t) + ^(*), (18) 

where £ a ,j/>W i s a rea ^ Gaussian white noise of correlation (£ a .i/i(^)£a,i/>(i')) = 
2D a 5(t — t') (same variance as £ a (t))- We can use Eq. (18) to obtain the Fourier 
spectrum ^(w) of the phase tp(t), and then its power density spectrum following 

(2Q)- 1 fj m (u) + \A0\- 1 LA") 2 



|*(w) 



H - 



(19) 



Note that here, the influence of gain on phase noise spectrum is not explicit 
anymore: it is implicitly contained in \Aq\. Figure 5 displays the phase noise 
spectrum explicitly expressed by Eq. (19), and we can now analyze how the 
spectrum behaves according to the various frequency ranges. 



5.1 Phase noise close to the carrier (u> < FIh) 

Here, we consider the spectrum for frequencies which arc relatively close to 
the carrier (with w > Ql, however). Qualitatively, this corresponds to the 
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Figure 4: An experimental power density spectra of the microwave noise below 
threshold. The spectrum has been scaled to its maximum (reference at ~ dB). 
(a) Spectrum in a 200 MHz window, showing how the noisy power density is 
profiled by the RF filter, (b) Zoom-in with near the 10 GHz central frequency 
in a 500 kHz window, showing the noisy ring-cavity peaks. 



frequencies that are much smaller than the high corner value Q.h of the multi- 
plicative flicker noise. In this region, flicker noise is stronger than white noise, 
so that \f] m {uu)\/2Q ^> \£ a ^(u))\. On the other hand, we can also consider that 
1 -e~ tuT ~ iuT. Therefore, taking into account the fact that [iT > 1, Eq. (19) 
can be simplified into 



I* 



close 



HI 5 



\Vm(u)\ 2 



n H D r , 



4Q 2 T 2 uj 2 2Q 2 T 2 uj 3 ■ 



(20) 



Some remarks can be made at this stage. First, The Lesson effect is here very 
explicit: the phase noise spectrum decreases as f~ 3 due to the / _1 flicker noise 
[9]. Secondly, The phase noise is inversely proportional to Q 2 . Thirdly, the 
phase noise is practically independent of the microwave amplitude \Aq\, as long 
as multiplicative noise is stronger than additive noise near the carrier. Finally, 
the phase noise decreases as T~ 2 , therefore justifying the need for very long 
delay-lines to reduce phase noise close to the carrier. This dependence was also 
recovered analytically by Yao and Malcki, using another theoretical approach 
[1]. Hence, in first approximation there are three ways to reduce phase noise 
close to the carrier: reduce the power D m of the flicker noise, increase the delay 
T or increase the Q factor of the RF filter. 
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5.2 Phase noise in the spurious peaks range (fl H < u> <C 
/*) 

The frequencies of concern are here those which are well within the bandwidth, 
but not too close from the carrier. This range typically lies between 50 kHz and 
few MHz, and contains the parasite ring-cavity peaks. It is also an area where 
the multiplicative and additive noises are both white [in the sense that in that 
range, |?7 m (cj)| 2 and l^a.v^)! 2 are both constant]. 

The local minima of the spectrum in that area are obtained for e~ lujT = — 1, 
so that the floor of the phase noise after the flicker decrease is 



I* floor? - \ 



2<9 \a | 



n 2 



(2D 



when Q is sufficiently high. This level is 6 dB below the additive white noise 
power density scaled to the power of the microwave. The recipe for a low phase 
noise floor is then quite simple, and also quite conventional: low additive noise 
D a , and high power |*4o| 2 for the microwave signal. The bandwidth does not 
play any role in this case, as long as the multiplicative noise is not too strong. 

The spurious peaks are localized around integer multiples of the round-trip 
frequency Q.t/2-k = l/T = 50 kHz. More precisely, a fourth order Taylor 
expansion of the denominator of Eq. (19) shows that around these resonance 
frequencies, the phase noise can be expressed as 



|#(nfi T + 5u)\ 2 = (22) 

2 

{2Q)- 1 f jrn {u) + \Aq\~ 1 Za,1>(">) 



M 2 



By finding the minima of this Taylor-expanded denominator, it can be shown 
that the spurious peaks are in fact frequency-shifted according to 

f n = - = nx50kHz-nxl6Hz. (23) 

T nAFT 2 

Then, their height relatively to the phase noise floor can also be calculated as 



A|ig^ B = 10 bg 



AF Tl 4 



n 



120 dB -401ogn. (24) 



It appears that the level of the spurious peaks increases with the RF bandwidth 
and with the delay: therefore, a large delay may lead to a lower phase noise 
near the carrier [see Eq. (20)], but it also leads to a higher level for the spurious 
peaks, so that an optimal trade-off has to be found. It is also noteworthy that 
this level is independent of the power densities D a or D m . In our case, the 
height of the first spurious peak relatively to the floor is theoretically equal to 
120 dB, in excellent agreement with the experimental results of Fig. 6, where a 
height of 119.5 dB has been measured. It can also be shown from Eq. (22) that 
the -3 dB bandwidth of the spurious peaks is 



(25) 
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Figure 5: A theoretical phase noise spectrum above threshold in a 500 MHz win- 
dow, with D a = 9.8 x 10~ 16 rad 2 /Hz, D m = 5x 10" 11 rad 2 /Hz, n H = 100 kHz 
and Q,l — 1 Hz. The dimcnsionlcss amplitude of the microwave oscillation is 
\Ao\ = 0.41, corresponding to a power of 10.5 dBm. 



an extremely small value which is experimentally confirmed with the results of 
Fig. 6b. In comparison, these spurious peaks typically have a linewidth higher 
than 1 kHz below threshold (see Fig. 3), but their linewidths sharply narrow as 
7-1. 



5.3 Phase noise outside the bandwidth (lj > fi) 

Here, the term /n[l — e~ luT ] progressively becomes negligible as u> is increasing, 
so that the ring-cavity peaks excited by white noise become strongly damped 
(for being outside the RF bandwidth). In this case, the phase noise decays as 



I* (c)\ 2 ~ 2Da ^ 1 ~ AF 2 ° a 1 (26) 

_ (26) 



However, the phase noise does not decrease monotonically as f~ 2 up to infinity: 
in fact, for lo 3> (J,, there is a second phase noise floor induced by the coupling 
between phase fluctuations and amplitude fluctuations (second-order effect, see 
ref. [4]). 



6 Conclusion 

This article has presented a theoretical study of phase noise in OEOs. Our ap- 
proach has consisted in a Langevin formalism, that is, in adding noise sources 
to a core deterministic model for the microwave dynamics. We have found a 
excellent agreement between the main predictions of the model and the exper- 
imental results. There is also an good agreement between this theory and the 
results that are known from the literature, or from our earlier works. 

The main advantage of this approach is that it enables within the same 
framework to understand the behavior of the system under and above thresh- 
old, as the same model continuously accounts for all the observed features in- 
dependently of the value of the gain. However, we have not taken into account 
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Figure 6: a) Experimental phase noise spectrum in a 100 kHz window, showing 
a noise floor around —145 dBrad 2 /Hz for a microwave power of P = 10.5 dBm. 
b) Enlargement of the spectrum around the first spurious peak at the frequency 
fi = 50594.35 kHz. The maximum of this spurious peak is at —25.5 dBrad 2 /Hz 
[height of the peak: —119.5 dB], and its —3 dB bandwidth is around 40 mHz. 
All these experimental data are in excellent agreement with the theory. Note 
that the height of the peak in Fig. a) is not indicative because of insufficient 
resolution. Also note that the peak at 50 Hz is a parasite peak originating from 
the electric mains supply. 



in this first model the noise generated by the filter (noisy /i and z9), and the 
delay time (noisy T). Fluctuations associated to these parameters may induce 
interesting stochastic features, that will be adressed in future work. Another 
line of investigation is to achieve a better spectral and statistical fitting of the 
multiplicative noise rj m (t), which is an essential variable for the determination 
of phase noise spectra. Future work will also emphasize on phase noise reduc- 
tion methods, such as optical filtering, multiple-loop architectures, or quadratic 
crossed nonlinearities [11]. 

A Determination of the output noise power for 
7 = 

In the open-loop configuration, the total output power can also be obtained 
using some quantum electronics formulas. 
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Effectively, the output power can be explicitly expressed as 

Po = \ [FkT + 2eI ph R eq ] GAF , (27) 

where G is the total gain of our two cascaded amplifiers (22.5 and 22.3 dB at 
10 GHz), F = 6 dB is the noise figure of the first amplifier, T = 295 K is 
the room temperature, k is the Boltzmann constant, e is the electron charge, 
I p h = 1.2 mA is the photodiode current, R eq — 25 ft is the equivalent load 
impedance for the photodiode, and AF = 50 MHz is the bandwidth of the RF 
filter. The formula gives Po = 19.6 nW, while we have measured 20.0 nW. The 
combination of Eqs. (16) and (27) also gives a method to determine D a directly 
from the specifications of the various optoelectronic components used in the 
oscillation loop. 



B Derivation of the stochastic phase equation 

We use ltd chain rules to derive the stochastic differential equation for the phase. 
We first rewrite Eq. (17) under the differential form 

dA = -iie^Adt + pe^[l + Vm (t)}A T dt + ^e vd dW a , (28) 

where dW a (t) — dW ar (t) + idW a i(t) is a differential Wiener process. Note 
that (dW a , r ) = (dW ai i) = and ((dW , r ) 2 ) = ({dW a ^) 2 ) = 2D a dt . The fact 
that dW ~ O(Vdt) explains why the differential terms of second order should 
be taken into account in stochastic calculus, so that usual differentiation and 
chain rules do not generally apply. When considering the second order one may 
consider (dW a ) 2 ~ ((dW a ) 2 ) = 2D a dt, and discard higher order terms since 
((dW a ) k+2 ) m 0[(dt) k ] < dt for k > 0. 

We set A(t) = ^fpjtj e^W = e p(t)+#(*) j where p= \ In V is an auxiliary 
variable. At order dt, we have 



,i t d -A 1 

dp + idip = d In A = — — - 



dA 



-i 2 



A 



+ ^ r dW a , (29) 

where we have considered {dW a ) 2 ~ ((dW a ) 2 ) = 0. Assuming second order 
fluctuations for the amplitude (that is, \At\ — \A\), we arc led to 

dtp = —fi sin $dt + p,[l + r] m (t)] sm{-d + tpT — ip] dt 

+ ^-dW a ^, (30) 

where dW a ^ = dW a , r sin $ + dW a ,i cos $ is a real Gaussian white noise with 
zero mean and variance ((dW a ^) 2 ) = 2D a dt. Since \ip — tpxl <S <S 1, we have 
sin[^ + ip T — tp] ~ i? — (ip — ip T ) so that finally, 

^ = -MV>-Vr) + ^?? ro (i) + ^W(*)- (31) 
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Figure 7: Oscillator phase noise transfer function (with s = jf). 



This result is also the one we may have recovered through the usual rules of 
differential calculus (however, note that it is not so for the equation ruling the 
power variable V). Also note that this equation is valid only as long as the 
approximation of neglecting B in Eq. (10) is valid. 



C An alternative paradigm for phase noise anal- 
ysis 

It is possible to gain a different physical insight into the phase noise problem in 
OEOs, using an alternative methodology related to the conventional theory of 
feedback oscillators. We hereafter briefly sketch the main lines of this heuristical 
approach. 

The oscillator consists of an amplifier of gain A (constant) and of a feedback 
path of transfer function (3{jf) in closed loop. The function (3(jf) selects the 
oscillation frequency, while the gain A compensates for the feedback loss. This 
general model is independent of the nature of the amplifier and of the frequency 
selector. We assume that the Barkhauscn condition \A/3(jf)\ = 1 for station- 
ary oscillation is verified at the carrier frequency / by through a gain-control 
mechanism. Under this hypothesis, the phase noise is modeled by the scheme 
shown in Fig. 7, in which all signals are the phases of the oscillator loop [9]. 
The main reason for describing the oscillator in this way is that we get rid of 
the non-linearity, pushing it in the loop-gain stabilization. The ideal amplifier 
"repeats" the phase of the input, for it has a gain of 1 (exact) in the phase-noise 
model. The real amplifier introduces the random phase ip(t) ^(jf) m the 
loop. In this representation, the phase noise is always additive noise, regardless 
of the physical mechanism involved. This eliminates the mathematical difficulty 
inherent in the parametric nature of flicker noise and of the noise originated from 
the environment fluctuations. 

The feedback path is described by the transfer function B(jf) of the phase 
perturbation. In the case of the delay-line oscillator, the feedback path is a delay 
line of delay T followed by a selector filter. The latter is necessary, otherwise the 
oscillator would oscillate at any frequency multiple of 1/T, with no preference. 
Implementing the selector as a bandpass filter (a resonator) of group delay T g , 
the phase-perturbation response of the feedback path is 

bo/) = c TlTf f T T) ■ < 32 > 



We assume that all the phase perturbations in the loop are collected in the 
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random function Tp(t) <-» ^f(jf), regardless of the physical origin (amplifier, 
photodetector, optical fiber, etc.). Denoting with <p(t) <-> ${jf) the oscilla- 
tor output phase, the oscillator is described by the phase-perturbation transfer 
function H(j/) = ®{j f) /^(j /)• By inspection on Fig. 7, and using the basic 
equations of feedback, the oscillator transfer function reads 

H(i/) = —L^, (33) 

and the oscillator phase noise spectrum would be given by S v (f) = |H(j/)| 2 S^(/). 
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